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In compensated elemental metals at low temperature, a several Tesla field can sup- 
press electronic heat conduction so thoroughly that heat is effectively earned by phonons 
alone. In approximately one mm diameter single crystal samples with impurity concentra- 
tions low enough that electron conduction is limited by surface scattering, the ratio of zero- 
field to high-field thermal conductivity can exceed ten thousand. We have used this phe- 
nomenon to build a compact, solid-state heat switch with no moving parts and no enclosed 
fluids. The time scale for switching states is limited by time scale for charging the magnet 
that supplies the controlling field. Our design and fabrication techniques overcome the dif- 
ficulties associated with manufacturing and assembling parts from single crystal tungsten. 
A clear disadvantage of the magnetoresistive switch is the mass and complexity of the 
magnet system for the controlling field. We have discovered a technique of minimizing 
this mass and complexity, applicable to the continuous adiabatic demagnetization refrigera- 
tor. 


INTRODUCTION 

Heat is carried in metals by both electrons and lattice vibrations (phonons). At low 
temperatures, the electron thermal conductivity is linear with temperature, while phonon 
thermal conductivity drops off with a higher power, and is typically several orders of mag- 
nitude smaller at liquid helium temperatures. A magnetic field suppresses the electron 
conduction in the plane perpendicular to the field. In metals that are compensated (those 
having an equal number of holes and electrons) and that have closed Fermi surfaces, the 
thermal conductivity drops with field until phonon conduction dominates. Metals for 
which this is true include Ga, Cd, Be, Zn, Mo, and W. Gorter and Miedema [1] realized 
that this effect could be used to make a low temperature heat switch, and implemented such 
a switch [2] with single-crystal gallium. In addition to its inconveniently low melting 
point, Ga has a rather high superconducting transition temperature, T c . For temperatures 
below T c thermal conductivity in the on (zero field) state drops as e c , severely restrict- 
ing the on-off ratio of the switch. If we are interested in a switch that operates well into the 
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subkelvin range, the choice of materials is quite limited. Radebaugh [3] pointed out that Be 
is perhaps the best material for a switch not only because of its extremely low T c (26 mK), 
but also because it has the highest DeBye temperature, @d, of any metal (1000 K). A high 
® D implies a very low phonon thermal conductivity, and indeed Radebaugh’ s measure- 
ments indicate that a Be switch would have a useful (> 10 3 ) on-off ratio up to 30 K. Unfor- 
tunately, it is very difficult to find large single Be crystals. Of the remaining candidates, 
tungsten stands out because of its low T c (15 mK) and reasonably high ©d (3 1 0 K). 

An important application for low temperature heat switches is in adiabatic demagneti- 
zation refrigerators (ADRs). ADRs have a number of advantages that make them ideal 
coolers for space applications, including high efficiency and independence from gravity 
(which allows full ground testing.) However, in a standard ADR a magnetoresistive heat 
switch would require a rather large magnet, to put it into the off state at the proper time. 
Because mass is such an important criterion for a spaceflight instrument, this requirement 
puts the magnetoresistive switch at a severe disadvantage relative to alternatives. 

The standard ADR has a problem of its own: to spend a large fraction of its cycle ab- 
sorbing heat at the operating temperature, the salt pill must be rather large. Simple scaling 
of current designs, such as the XRS ADR [4], to meet the requirements of future missions 
leads to systems with unreasonable size and mass. Recently, Shirron, et al. [5] have cir- 
cumvented this limitation with the invention of the continuous ADR (CADR). A CADR is 
a multistage ADR with the stages connected in series. Stage 1, which is attached to the 
load, acts as a normal ADR during demagnetization; the field is ramped down just fast 
enough to absorb the heat of the load. However, before this stage reaches zero field, stage 
2 is demagnetized to a temperature below the operating temperature and a heat switch be- 
tween them is closed. At this point, the stage 1 must magnetize to maintain the operating 
temperature. Once the stage 1 reaches full field, the heat switch is open and it starts de- 
magnetizing again. When the heat switch opens, stage 2 is ramped up. Simultaneously, 
stage 3 is ramped down. When the temperature of stage 2 rises above that of stage 3, the 
heat switch between them is closed, and stage 2 transfers its heat to stage 3. This cascading 
of heat from one stage to the next may be repeated by adding additional stages to raise the 
final heat sink temperature. The key point is that by decoupling the operation of the device 
being cooled (typically detectors) from the operation of die ADR, the ADR may be de- 
signed to cycle much more rapidly, thus providing higher cooling power unit mass. 

Note that for any pair of stages in the CADR, the magnetic field of the higher stage is 
low when the switch needs to be closed, and it is high when the switch needs to be open. 
Suppose the part of the thermal link between the two pills that is exposed to the magnetic 
field of the higher stage magnet is a magnetoresistive material. Then it appears that the 
material will act as a “passive” switch, opening and closing at the proper phase in the cycle 
without a separate control circuit. 


MAGNETORESISTIVITY OF TUNGSTEN 


Tungsten has a simple body centered cubic crystal structure. Because of this, if the 
field is applied along a high symmetry direction (e.g. 001), the thermal conductivity tensor 
will have the form 
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Measurements of the transverse conductivity X^ by Long [6] and Fletcher [7] indicate that 
it is never more then a few percent of X^ over the range of temperatures and fields of inter- 
est. Thus, thermal conductivity is isotropic in the plane normal to the field to within this 
error. 

These and other authors [8], [9] have investigated the electrical and thermal resistivity 
of single crystal tungsten at helium temperatures and fields up to several Tesla. All find 
that the component has the form 

(?) 


where Tg is the phonon conduction term. Wagner [8] found that in high field, A xx fits to 
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Fitting to his data for a very high quality (RRR = 63,000) crystal gives cto = 
0.42 Tesla 2 W m' 1 K 2 and a 3 = 0.0063 Tesla 2 W m' 1 K' 5 . The temperature independent 
part, due to impurity and surface scattering, is related to the temperature independent part 
of the electrical resistivity, by ao * To B^/pxxO, where Lq is the Lorenz number. 

Long [10] and Wagner [8] find that Tg is proportional to T 2 with constants of propor- 
tionality 0.038 and 0.050 W m 1 K 1 , respectively. This result is unusual in two ways: the 
dependence is quadratic, rather than cubic, as is typically found in dielectrics, and the val- 
ues in tungsten crystals are far lower than one would expect in a good dielectric crystal. 
Long shows that this is because phonons are scattering primarily off conduction electrons. 

In zero field, the thermal resistivity, w = k~ l , at low temperature is given by 

wT=(wT) T =o + a1 2 , (4) 

which has the same form as electrical resistivity, p(T) = po + AT 2 . Again, the constant 
term, {wT)t=o! po ~ To, as one would expect from impurity and surface scattering. On the 
other hand, Wagner et al. [11] find that a JA is much smaller than L 0 . For the same sample 
used to derive <xo and a 3 , (wT)t=o = 0.00316 K 2 cm/W and a = 2.21 x 1CU 4 cm/W. 

Unfortunately, there is no good theory explaining the low field region. For the pur- 
poses of modeling the behavior throughout the range of fields, we assume the magneto- 
thermal resistance acts in parallel with the zero-field resistance: 


k(T,B) = k(T)+- 


Using the above coefficient values in this equation, one finds that zero field to high 
field thermal conductivity ratios of over 10 5 are possible for T < 2K and B > 2Tesla. Un- 
fortunately, for heat switches in the CADR, high thermal conductance is needed at low 
temperatures, and in the “ofF’ (low conductance) state the switch must span a large tem- 
perature range, in which case the conductance is dominated by the thermal conductivity of 
the high temperature end. The important quantity, the ratio of the heat flows in the two 
states, is given by the ratio of the integrated thermal conductivity. At high fields (although 
not so high that k « k g ) and low temperatures, 
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Thus, at low temperatures, the ratio of heat flows is given approximately by 

Qon w 2 B T m , g on AT^ 

Qoff ( W7, )r=0 a 0 Thigh, off Thigh, off 
This has the advantage that the parameters ( wT)r=o ar >d «o can be estimated from electrical 
resistivity measurements. For the values given above, 7 = 3.65 mTesla. For a 

typical case: Thigh, on = 0-25 K, Ti ow , on = 0.225 K, Thigh,off = 1-3 K, Ti 0W , 0 jf = 0.045 K, B = 
1.0 Tesla, this equation gives a ratio of 472, whereas integration of Eq (5) gives 529. Inte- 
gration must be used at higher temperature; for the case Thi g h.on = 1-3 K, Ti ow , on = 1.2 K, 
T high ,off= 4.2 K, T low .off= 0.225 K, B = 2.0 Tesla, it yields a ratio of 1 .39 x 1 0 3 . 

For the CADR, heat leak through the switch in the off state, Q off , enters the colder salt 


pill at T i ow , 0 ff and must be rejected by heat flow in the on state, , at Thigh, on- Assuming 

that the pill spends half the cycle period, t, at Ti ow , 0 ff and half at Thigh.on, and that the cycle 
efficiency is 1 00% of Carnot, to do useful cooling, we must have 
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In fact, the left side, the enthalpy rejected in the on state due to heat leak in the off state, 
should be a small fraction of the right side, the total enthalpy rejected through the switch in 
the on state. If we demand that it be no more than 1 0%, then in the high and low tempera- 
ture examples above we find requirements of Q on / Q off > 56 and 58, respectively. Thus, 
very high purity tungsten easily meets this requirement in the two example cases. 


SWITCH DESIGN 

The area-to-length ratio of the switch, All, is set by the requirement that all the heat 
absorbed in the lower stage when it is at Ti ow o ff be rejected through the switch. Again as- 
suming a 50% duty cycle and 1 00% efficiency, this condition is the same as Eq. (9) with 
Qoff replaced by the total heat load on the lower stage. The on state heat flow is just 

All ( Ushm k(T,0)dT . Using a load of 25 pW and the temperatures in the previous example, 

*lvw,on 

and allowing a factor of two margin to account for inefficiencies gives All = 5.2 x 1 0 m. 

In experiments on samples of very high purity (RRR « 75000), Batdalov [9] found the 
electron mean free path, 1.4 mm T = 4 K. Reducing the dimensions of the conductor 
below d significantly increases the effective resistivity. Further, making the width and 
thickness of the conductor very small also makes the structure overly fragile. For these 
reasons, a square cross section of 1 .5 mm was chosen. This implies that an effective length 
of 0.43 m is required. 

As mentioned in the introduction, the key to making the magnetoresistive switch at- 
tractive for use in a CADR is that it can exploit the “existing” field to turn the switch off. 
However, this is not quite true. To accommodate the switch, the magnet must be longer or, 
equivalently, the salt pill must be smaller. This runs counter to the prime advantage of the 
CADR: its high cooling power per unit mass. Thus, we need to minimize the length of 
magnet bore that the switch occupies. 

Our solution is shown in FIG 1. Starting with the largest available tungsten single 
crystal, we wire EDM (electronic discharge machine) a disk 20 mm in diameter and 5.1 
mm thick. The cylinder (z) axis is aligned with the 001 direction. The other crystal direc- 
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FIGURE 1 . Design of magnetoresistive heat switch. The shape is wire EDM cut from a tungsten single 
crystal. 

tions were not known, but, as discussed above, conductivity in the xy plane should be 
nearly isotropic. The stated purity was 99.995% . Two slots are then cut in the horizontal 
(xy) plane, and 9 in the xz plane, as shown. The result is effectively a 1.5 mm x 1.5 mm 
rod, approximately 0.4 m long, including only the length of the straight regions. A flexible 
thermal strap of high purity copper wires is vacuum brazed to the “D” shaped end section. 


MEASUREMENTS 

We used x-ray diffraction to characterize the structure of the crystal. A small sample 
taken from a wire EDM machining test was measured. Nearly all of the power appeared in 
a single Bragg peak, indicating the material was indeed a single crystal. The approximately 
one percent of power that appeared randomly aligned was attributed to small surface ciys- 
tallites caused by rapid remelting during wire EDM machining. 

To measure the magnetoresistance of the material, a piece similar to the one shown in 
FIG.l, but with only a single layer, was cut from the crystal. Four copper leads were vac- 
uum brazed to the sample. At room temperature, the resistance of the sample was 0.0035 
Ohm. Using Long’s [10] value for room temperature resistivity of pure tungsten, 5.5 mi- 
croOhm-cm, gives All =1.61 0" 5 m, which is 7% larger than that calculated, considering 
only the straight sections of the sample. The sample was mounted on a probe that fit in the 
bore of a H^p«*Supervarimag cryostat. A 0. 1 Amp peak-to-peak ac current was applied 
and the voltage was read on with a lock-in amplifier. Measurements were made at 7'= 4.26 
K. The resulting resistivity as a function of field is shown in FIG 2. The zero-field resis- 
tance was so low that parasitic reactances dominated the signal into the lock-in. Using a 2 
Amp dc current and a sensitive voltmeter, we were able to place a lower bound on the zero 
field resistance of 0.5 microOhm, implying a lower bound of 7 x 10 3 on the RRR. 

Thermal conductance of the material was measured as a function of field using an ear- 
lier version of the switch. The geometry was the same as shown in FIG. 1, but the bottom 
layer was much thicker than the two “active” layers. A ring of commercially pure 
(99.95%) tungsten was brazed around the bottom layer creating an approximately isother- 
mal surface, to which a ruthenium oxide thermometer was varnished. Two bare ruthenium 
oxide chips were varnished onto the slender arm of the other side of the switch. One was 
used as a heater. The other was calibrated in situ and used as a thermometer. The switch 
was placed in a vacuum can that fit into the Mtmm Supervarimag. A heavy copper braid 
was used to link the bottom side of the switch to the flange of the vacuum can. 





=1 



FIGURE 2. Resistivity as a function of field. The solid line is a power law fit. 

The thermal resistivity, w, times temperature verses magnetic field is plotted in FIG 3 
for several different temperatures. The data fit well to a power law: wT = wT\Ep with p = 
1.7. The fact that the exponent is not 2 is somewhat surprising. However, the electrical 
resistivity verses magnetic field data also fits to a power law with an exponent of 1.7. The 
two samples have substantially different geometries. Long points out that his data actually 
fits best to an exponent of 1 .96. Fitting the coefficients wT\ to the form of EQ.(3) yields a 0 
= 1.37 W Tesla 1 7 m' 1 K‘ 2 and a 3 = 0.0061 Tesla 2 W m* 1 K 5 . This value of a 3 is essentially 
the same as obtained from Wagner’s data. The value of oto is larger, as one might expect 
from material with higher impurity content. Fitting the zero thermal resistance verses tem- 
perature to the form of EQ.(4) yields (wT)jh>= 0.029 K 2 cm/W and a = 2.3 x KT 4 cm/W. 

These coefficients can be used to compute Q m / Q off for this material. Substituting the 

coefficients into EQ.(5), changing the exponent on B to the measured value, and integrating 
yields ratios of 17 and 78 for the high and low temperature cases discussed previously. In 
the low temperature case, the ratio does not meet the requirement. A field >2 Tesla, higher 
than desirable for this temperature range, or a higher purity material is needed. However, 
the ratio in the high temperature case exceeds the requirement (58) by a significant margin. 


MODELING SYSTEM PERFORMANCE 

A simple ADR model was modified to simulate a two stage ADR with a “passive” 
magnetoresistive switch. That is, the field applied to the switch is assumed to be the same 
as that applied to the higher temperature (“upper”) salt pill, which is also connected to a 
fixed-temperature bath through a gas-gap heat switch. The on conductance of the gas gap 
switch is realistically modeled, but off conductance is ignored. The stages below the lower 
stage are simulated with a fixed heat load if the pill temperature was near its lower setpoint, 
and zero otherwise. The entropy and heat capacity of the salt are modeled using the stan- 
dard equations for simple, non-interacting magnetic ions. Interactions are taken into ac- 
count by replacing the applied field, B, with V B 1 + b 7 [12]. The “interaction field , b, for 
various paramagnetic salts is derived from fitting to zero field entropy data. Despite the 
simplicity of the entropy function, this model works quite well for single stage ADR’s, par- 
ticularly when using salts with small crystal field splitting and low Neel temperature such 
as FAA (Ferric Ammonium Alum) and CPA (Chromium Potassium Alum). 





FIGURE 3. Thermal resistivity times temperature of the tungsten sample as a function of magnetic field at 
several temperatures: x, 1.6 K; +, 2.2 K; □, 4.0 K; O, 8.0 K. The line is a power law fit to the 2.2 K data. 

Simple proportional controllers set the voltage across the magnets. The error signal is 
to the difference between the temperature and a setpoint. For the lower stage, the algorithm 
for deciding the setpoint is simple: When the stage is in the magnetizing state, a “high 
temperature setpoint is used. Once the magnetic field approaches the maximum allowed 
value, the stage switches to the demagnetizing state with a “low” temperature setpoint. As 
the field approaches zero, the stage returns to the magnetizing state. The setpoint algorithm 
for the upper stage is the same for the magnetizing state. In the demagnetizing state, how- 
ever, it is somewhat more complex. If the lower stage is still demagnetizing, the setpoint is 
put at an intermediate value. Although the field on the switch is lower, the temperature dif- 
ference is also, and a value can be chosen which minimizes heat leak. Once the lower stage 
is ready to magnetize, the setpoint temperature is calculated to give a fixed heat flow across 
the switch, unless a minimum temperature is reached. 

FIG 4 shows the results of one simulation. In this run, the upper and lower pills were 
FAA and CPA pills with masses of 1 50 and 50 gm, respectively. The plot shows the third 
cycle in the simulation, after initial transients have settled out. The first part of the plot is 
uneventful. The lower stage temperature rapidly reaches its lower setpoint, and the field 
gradually decays due to the applied heat load. Once the lower stage starts magnetizing, a 
significant issue for the passive magnetoresistive switch becomes clear: To achieve the de- 
sired heat flow in the on state, the upper salt pill must be demagnetized to a temperature far 
below the upper setpoint of the lower stage. Absorbing the heat at lower temperatures re- 
sults in lower thermodynamic efficiency. For the cycle shown, the calculated efficiency, as 
a fraction of Carnot efficiency, is ( QioaJQbath ) x ( Tbath/Ttoad ) = 32%. Approximately one 
quarter of the efficiency loss is due to heat leak through the magnetoresistive switch in the 
off state. This can be greatly reduced by operating with a higher maximum field in the up- 
per stage, but at a loss in overall thermodynamic efficiency because the salt pill must drop 
to a lower temperature. 

One way to greatly reduce this loss is to choose a paramagnetic material for the upper 
stage with an ordering temperature closer to the maximum temperature of the lower stage. 
In this case, the temperature of the upper stage cannot drop far below that of the lower 
stage, and the field drops rapidly as the temperature approaches the ordering temperature, 
allowing the system to take full advantage of the switch’s high low field conductance. 





FIGURE 4. Calculated temperature and field behavior for a two stage ADR with a magnetoresitive switch 
between the stages. 

CONCLUSION 

Heat switches based on the thermomagnetoresi stance of high purity metals have po- 
tential application in ADRs. They are particularly attractive for the continuous ADR be- 
cause the field of each stage could control the thermal resistance of the switch to the lower 
stage, thus eliminating the cost and complexity of a separate control magnet. The perform- 
ance of such a switch, calculated using material properties from the literature, indicate that 
it should easily satisfy the requirements of the CADR. A switch of the appropriate geomet- 
ric properties was wire EDM machined from a large tungsten crystals. Measurements on 
samples cut from this crystal indicate that it will have somewhat lower, although still ac- 
ceptable, performance. A model was developed to simulate a two-stage ADR using a mag- 
netoresistive switch. Though relatively simple, the model points out that true system per- 
formance depends in a complex way on the interactions between the fields, temperatures, 
and the control algorithm. It shows that if the ordering temperature of the salt in the higher 
temperature pill is too low, the system will lose efficiency because the pill will need to drop 
to a lower than desirable temperature to achieve sufficient interstage heat transfer. 
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